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The random-phase-approximation semiclassical scheme for description of plasmon excitations in 
large metallic nanospheres, with radius range 10-60 nm, is formulated in an all-analytical version. 
The spectrum of plasmons is determined including both surface and volume type excitations and 
their mutual connections. The various channels for damping of surface plasmons are evaluated and 
the relevant resonance shifts are compared with the experimental data for metallic nanoparticles of 
different size located in dielectric medium or on the semiconductor substrate. The strong enhance- 
ment of energy transfer from the surface plasmon oscillations to the substrate semiconductor is 
explained in the regime of a near-field coupling in agreement with recent experimental observations 
for metallically nanomodified photo-diode systems. 

I. INTRODUCTION 

Rapid progress in plasmonics 1 (taking advantage of peculiar properties of plasmon-polaritons^ in nano-structured 
metallic interfaces) and plasmonic applications in photonics and microelectronics^, focused attention on metallically 
modified systems in nano scale and on collective excitations of metallic plasma in a confined geometry. Of particular 
interest are recently reported experimental data on giant enhancement of photoluminescence and absorption of light 
by semiconductor surface (of photo-diode) covered with metallic (gold, silver, or copper) nanospheres, with sphere 
radius of order of several to several tens of nanometers^&£&2. 

These phenomena are considered as perspective for enhancement of efficiency of solar cells by application of special 
metallic nanoparticle coverings of photo-active layer^. Metallic nanospheres (or nanoparticles of other shape) can 
act as light converters, collecting energy of incident photons in surface plasmon oscillations. This energy can be 
next transferred to semiconductor substrate in a more efficient manner in comparison to the direct photo-effect. 
The experimental observations^££>£>£>2 suggest, that the short range coupling between plasmons in nanospheres and 
electrons in the semiconductor substrate allows for significant growth of selective light energy transformation into 
a photo-current in diode system. This phenomenon is not described in detail as yet, moreover some competitive 
mechanisms apparently contribute, manifesting themselves in strong sensitivity of the effect to size and shape of 
metallic nanocomponents, type of the material and dielectric coverings of nanoparticle o 9 ' 10 . It is probably connected 
with particularities of the short range dipole-type (similar to the Forster interaction}^ coupling of sphere surface 
electrical oscillations with semiconductor substrate. Nevertheless, one can argue generally that due to nanoscale of 
the metallic components the momentum is not conserved, which leads to the allowance of all indirect optical intcrband 
transitions in semiconductor layer, resulting in enhancement of a photo-current in comparison to the ordinary photo- 
effect when only direct interband transitions were admitted. Though the effect is selectively ranged to a vicinity of 
resonant plasmon frequency, the gain in the total efficiency would be enlarged by possible dispersing of dimension 
and shape (or by a dielectric coating) of metallic nanoparticles, widening the resonant spectrum. Any efforts towards 
improvement of the efficiency rate of photo-cells via not complicated metallic modifications of their surfaces are of 
particular significance, as the main drawback on the way for a wider application of these renewable electricity sources 
is their relatively small efficiency (reaching not more than 10-20%, depending on a material). 

Since in the metallically modified photo-cell structures the surface plasmons play a central role, the recognition 
of these excitations in nanoparticles is important. The surface plasmons have been originally considered by Mie^i, 
who provided a classical description of oscillations of electrical charge on the surface of the metallic sphere within the 
electron-gas model. The dipole-type Mie surface oscillations of the electron plasma are not dependent of the sphere 
radius, in contradiction to experimental observations both for small metallic clusters and for larger spheres. The 
plasmon effects are of particular significance in noble metals (gold, silver, and also copper) due to strong visible-light 
plasmon resonances in these materials. Even though a bulk luminescence rate in metals is small (~ 10~ 10 ), the 
rough surface structured in the nano-scale exhibits luminescence enhancement by several orders of magnitude. A 
rough metal surface can be treated as a collection of randomly oriented nano-hemi-spheroids, with distinct spectra 
of surface plasmons, increasing the luminescence rate (e.g., by 6 orders in golcU 2 -) — in particular it leads to an easily 
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observable 'multi-colored lightning rode effect— on the surface of broken polycrystalline metallic sample. Note also 
that triangle-branched recently synthesized gold colloidal nanoparticlesi 4 - are blue in contrast to spherical ones being 
red, which is caused by the shape-induced shift of surface plasmon resonance^ 5 -. 

In order to describe all these phenomena the plasmon properties beyond the Mie-type picture need to be recognized. 
Plasmon excitations in metallic clusters were the subject of wide analyse s 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 ' 24 within various attitudes 
including quantum effects. The quantum effects in small metallic clusters were analyzed using numerical methods of 
calculus 'ab initio', including shell-model and Kohn-Sham 'local density approximation' [LDA], similar as applied in 
chemistry for large molecule calculations, limited however to few hundreds of electron s 16 ' 18 ' 19 ' 20 ' 22 . Also variational 
methods for energy density and the random-phase-approximation (RPA) numerical summations were applied (e.g., for 
clusters of Na with radius ~ lnm) 2 ^ and various semiclassical expansion method s 17 ' 20 ' 21 . Emerging of the Mie response 
from the more general behavior was presented, in particular, decoupling of surface — translational and volume — 
compressional modes in small clusters was demonstrate d 16 ! 18 ' 25 . In description of metallic clusters commonly the 
'jellium' model was adopted, allowing for adiabatic approach to background ion system. In jellium model all the 
kinetics concerns electron liquid screened by static and uniform positive background of ion o 16 ' 26 ' 27 ' 28 . 

Problem of plasmon oscillations in metallic clusters was analyzed also including analytical formulation employing 
Thomas- Fermi- type approach, e.g., by Kresin 1 — or more recently within semiclassical expansion and separation of mass 
center and relative electron dynamics 2 ^, without, however, inclusion of damping of surface plasmon oscillations due 
to radiation losses and addressed rather to small systems (up to ca 200 atoms, i.e., at most 2.5 nm for radius 2 ^). The 
similar range of sphere size was assumed in more microscopic approach, time-dependent local density approximation 
(TDLDA), widely applied to metallic cluster analyses (e.g., Refs [T^)l8|f20h . where the scheme above Thomas-Fermi 
approximation was developed, including spill-out effect of electron cloud beyond positive jellium, screening effects 
and decay of surface plasmons into particle-hole pairs (Landau damping) 19 ! 21 . TDLDA type of approach, though 
also basing on jellium model, in a much better way allows for experimental data elucidation, especially for small 
and ultra-small clusters, with radius up to 1-2 nm, mostly of Na or K and noble metals (however, still not with 
all particularities) 16 i 18 . In this range of cluster dimensions the experimentally observed red-shift of Mie frequency is 
described mainly by spill-out effect, which is pronounced in small nanocrystals and causes a red-shift due to reducing 
of electron density (Mie frequency is proportional to the square of the density) 17 ! 18 ' 20 . But for dimensions above only 
few nanometers, the lowering of electron density caused by spill-out is of order of single percent or smaller—, and thus 
for higher nanosphere dimensions (of order of few tens nm) the spill-out effect is rather unimportant (the resonance 
frequency shift caused by spill-out is proportional to inverse radius of a metallic sphere, as being typical surface-type 
effect, and thus diminishes with radius growth). The other quantum effects, as magic electron numbers, related with 
closed shells for confined syste m 17 ! 18 correspond also to small clusters, starting from only several electrons to few 
hundreds of electrons (though supershell beats, up to N = 3000 in alkali clusters, were also investigated — for review 
cf. Ref. Il6h . For ultra-small clusters electron collective excitations are strongly coupled between volume and surface 
and also with background ionic system oscillations. Emerging of a well formed volume plasmon for ^0.5 nm radius 
was demonstrate d 16 ' 18 and it was argued that for smaller clusters separation of surface and volume excitations is 
impossible (calculations are ranged to 1=1 dipole mode 1 ^). With growth of the cluster radius the situation changes 
and both types of collective excitations: translational and compressional, addressed to surface and volume oscillations, 
respectively, can be separated (finally they have strongly different frequencies). Note, that useful notions of volume 
and surface oscillations were also applied in the case of nuclear matter vibrations of compressional typ o 29 ' 30 and 
of translational type 3 ^, respectively, which was then important in understanding of giant nuclear dipole resonance 
experiments. 

The Landau damping due to decay of plasmons for particle-hole pairs with similar as plasmon energy was also 
analyzed (in order to explain observed experimentally resonance frequency red-shift for small clusters) which is, 
however, also a diminishing effect with growing radiu s 19 ' 24 . Decay of surface plasmons into particle-hole pairs (Landau 
damping) leads to red-shift of Mie frequency (not monotonic with radius for ultra-small clusters)^ and similar in 
scaling, ~ i, to spill-out-caused shift, for larger sphere s 19 ' 21 (Landau damping is rather limited to radius range up to 
2.5 nm)2&. 

In the present paper we develop the RPA semiclassical method, originally formulated by Bohm and Pines for bulk 
meta l 32 ' 33 , in order to describe electron collective excitations in large metallic nanosphere, with radius range 10-60 
nm (much larger than Thomas- Fermi radius being of order of interparticle separation), including both volume and 
surface types of plasmons in the framework of an all-analytical calculus. The metal is assumed as a so-called 'simple 
metal', i.e. allowing for description of the electron-ion interaction by local and not strong pseudopotential (this 
condition is satisfied e.g., for noble, alkali, or transition metals) 32 . In the next section, the RPA equations for a local 
electron density are derived, including conditions imposed by finite geometry of the nanosphere (particularities of 
the solution method are shifted to Appendix [X| . In the following section spectra of surface and volume plasmons 
for the metallic nanosphere are presented, along with their modification by dielectric medium in which a metallic 
sphere can be embedded (the surface plasmon frequencies exhibit significant dependence on the dielectric constant 
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of a surrounding material in contrary to the volume plasmon frequencies). Within this semiclassical RPA attitude, 
the influence of the volume excitations onto the surface ones is also described, including possible exciting of surface 
plasmons by volume electron fluctuations. Finally, the e-m response of dielectric medium with metallic nanosphere 
subsystem is analyzed in the case of the dipole type (1 = 1) excitation, including modifications caused by plasmon 
damping, resulting in a strong dependence of resonance energy shift on a nanosphere radius, which was observed 
experimentally for large nanospheres^. Various channels of plasmon damping are considered including radiation 
losses in far- field and near-field regimes (in Appendices B and C, respectively) and the resulting resonant spectrum 
modification is compared with experimentally measured dependence of emission and absorption rates with respect to 
the sphere radius^ and dielectric coating^^. The giant strengthening of a coupling in a near-field regime of surface 
plasmons with semiconductor substrate is described in agreement with the experimental data for enhancement of a 
photo-current in metallically nanomodified diode system^jSiZiSiS. 



II. RPA SEMICLASSICAL APPROACH TO ELECTRON DISTRIBUTION IN METALLIC 

NANOSPHERE 

A. Derivation of RPA equation for local electron density in a confined geometry 



Let us consider a metallic sphere with a radius a located in vacuum, e = 1, p, = 1. We assume that the interaction 
between electrons and ions is described by a local and weak pseudopotential (this condition corresponds to the so-called 
'simple metal' case)22, as e.g. for noble metals. The Hamiltonian for this system has the form: 



^U 2 v 2 i x . ,„ 

H = - > v - + - > u(R v 
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where i2„, Vj and M, m are positions and masses of ions and electrons, respectively; N — number of ions in the 
sphere, N e = ZN — number of collective electrons, u(R u — R u ') is the interaction of ions (ion is treated as a nucleus 
with electron core of closed shells), w(R u ~ Vj) is the local pseudopotential of electron-ion interaction. Assuming 
the jellium mode l 16 i 26 ' 27 one can write for the background ion charge uniformly distributed over the sphere: n e (r) = 
n e <d(a — r) where n e — N e /V and n e \e\ — averaged positive charge density, V — ^p- — sphere volume, is the 
Heaviside step-function. Then, neglecting ion dynamics and small electron-ion pseudopotential (shifted by jellium- 
electron interaction), collective electrons can be described by the Hamiltonian: 
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with corresponding electron wave function \E' e (t). 

A local electron density can be written as follow o 32 ' 33 : 



P (r,t) =<* e (t)\y,t(r-r j )\* e (t)>, 



(2) 



(3) 



with the Fourier picture: 



where the 'operator' p(k) = ^ 



p(k,t) = / p(r,t)e- lkr d d r =< * e (t)|p(fc)|* e (t) >, 
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Using the above notation one can rewrite H e in the following form, in an analogy to the bulk casi 
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where: n e (k) = / d 3 rn e (r)e lk v , |j = J d 3 r^e 

Utilizing this form of the electron Hamiltonian one can write out: 
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in the following form: 



X>- ifc -' -v,) > 

If one takes into account that /3(fc — q)p(q) — 6p(k — q)Sp(q) + h e (k — q)8p(q) + 5p(k — q)n e (q) + n e (k — q)n e (q) 
and h e (q)p(k — q) = n e (q)5p(k — q) + h e (q)h e (k — q), where 5p(k) = p(k) — n e (k) describes the 'operator' of local 
electron density fluctuations above the uniform distribution, one can rewrite Eq. in the form: 

3 1 J (8) 

/ d 3 qn e (k " q)^Sp(q) - ^ J d?q5p{k - q)^6p(q). 

Thus for the electron density fluctuation, Sp(k,t) —< ^ e \Sp(k, t)\^> e >= p(k,t) — n e (k), we find: 

2^ = -<* e | Ee -^{-^ (fc .v,) 2 + ^fc.v, + ^}|* e > 

j (9) 
~M? S d3 lMk - q)^Sp(q, t) - ^ J cfq^- < ^> e \5p{k - q)8p{q)\* e > . 

Within scmiclassical approximation three components of the first term in right-hand-side of Eq. @ can be estimated 

as: k 2 Vp5p(k,t), k 3 vp/kTSp(k,t) and k i v 2 F /k 2 ^5p{k,t), respectively (1/fcr is Thomas-Fermi radius 32 , &t = ^/ 67r ^ e2 ; 

tp — Fermi energy, vp — Fermi velocity). The contributions of the second and the third components of the first term 
can be neglected in comparison to the first component. Small and thus negligible is also the third term in right-hand- 
side of Eq.®, as involving a product of two dp (which we assumed small, Sp/n e << 1). This approach corresponds 
to random-phase-approximation (RPA) attitude formulated for bulk meta l 32 i 33 (note that (5/5(0) = and the coherent 
RPA contribution of interaction is comprised by the last but one term in Eq. ([9"|)). 
Within the RPA Eq. © attains thus the form: 

d 2 Sp(k,t) 2k 2 lNr ^ Ti 2 V 2 „ Aire 2 f ~ , ^fc•<Z,,^ 

where for the case of spherical symmetry: 

h 2 9k 2 ?i 2 V 2 

< * e | V e~ ik ^ — (k ■ V,) 2 |* e >~ — < *«| V e~ ik - r * — -L\9 e > . 

m z 3m ' 2m 
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In the position representation Eq. (|10|) can be rewritten in the following manner: 



9I W 1 = < *-l E «(r - rj)^|*. > 

3 

+ gv {e(a - r)V / d 3 r 1F ^<5p(r l7 <)} . 
The Thomas-Fermi averaged kinetic energy can be represented as follows^ 2 -: 

<* e |-5: < J(r-r J )^|* e >^|( 37 r 2 ) 2 / 3 |i(p(r,t))5/ 3 



(11) 



3 r n (12) 

| (3 ^ ) 2/3£ n 5/3 e(a _ r) 
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Note, that neglected here gradient terms (in particular von Weizsacker term, ~ (Vp) 2 /(4p), beyond the Thomas-Fermi 
formula for kinetic energy functional ~ p 5 / 3 ^I& strongly affect the finite system properties especially of small metallic 
clusters. The contribution of this particular term (von Weizsacker) depends on the approximation in various versions 
of corrections to Thomas-Fermi approach 1 - (the coefficient of von Weizsacker term is treated even as a convenient 
fitting parameter). The gradient terms are inexplicitly included in TDLDA type methods based on Kohn-Sham 
equation. As it follows from respective analyses the contributions related to these terms (mostly spill-out effect) are 
more important for small clusters (when the surface dominates) and gradually diminish with the nanosphere radius 
growthi^^^^i. 
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Taking then into account that V0(a — r) = ——S(a — r) = — £ lim e _>o 5 (a + e ~ r )7 one can rewrite Eq. (fTTjl in the 
following manner: 

^§^= [|^V 2 <5p(r,t)-^p( r ,t)]e( a -r) 

-^V{[|e F n e + eF ^(r,f)] £% + e-r)} (13) 
!f^V% i) + fev / d 3 r 1F ^^(r l5 «)] <5(a + e - r). 

In the above formula lo p is the bulk plasmon frequency, lo^ = 47r ^ e , and the abbreviated notation, 5 (a + e — r) — 
lim £ ^ ( 5( a + € — r )i was used. The solution of Eq. (fT5|) can be decomposed into two parts related to the distinct 
domains: 

5Kr,t) = {¥f^ f f ° rr < a > . (14) 
|_ op 2 ir, t), for r > a, (r — > a+), v ' 

corresponding to the volume and surface excitations, respectively. These two parts of local electron density fluctuations 
satisfy the equations: 



and 



— = - 3^ V { [hpn e + e F 5p 2 {r, *)] ^8{a + e - r)} 



W 2 5p 1 {r,t)-J 1 p 5p 1 {r,t), (15) 



2 Sp 2 (r,t) _ 



^VSp 2 (r, t) + i^V / d 3 n ^ (Sprin, t)G(a - n) + <J/5 2 (n, t)0(n - a))] 5(a + e - r). 



Within this quasiclassical simplified approach volume plasmons described by the Eq. (|15p are independent of the 
surface plasmons, while the latter can be excited by the former ones, due to the last term in Eq. p6[) (which is 
caused by a 'surface tail' of compressional — volume type oscillations, while oppositely, translational — surface type 
oscillations do not have a 'volume tail'), which expresses a coupling between surface and volume plasmons in large 
metallic nanosphere. Coupling between volume and surface plasmons was analyzed in TDLDA approach for jcllium 
sphere s 16 ' 18 and it was demonstrated that for ultra-small clusters it is impossible to decouple volume and surface 
oscillations (because of perturbation of a shell structure), while for bigger clusters (more than 58 electrons, for Na 
cluster) the well formed and separated both modes emerg o 16 ' 18 (the TDLDA analyses were done for clusters up to ca 
200 electrons). This volume-surface coupling is strong for ultra-small radii, when shell effects and quantum spill-out 
are pronounced, and gradually weakens with growing sphere dimension. 

In the present paper we consider the radius range 10-60 nm, when quantum effects, significant for smaller clusters, 
are not of primary importance and are dominated by irradiation behavior. Such large nanosphcres contain 10 5 -10 7 
collective electrons, thus are considerably larger than clusters with up to 200 electrons, numerically investigated in 
details by use of TDLDA methods. Our main idea is to formulate analytical (thus simplified) RPA description in 
the form of oscillator equation allowing for phenomcnological inclusion of damping rates due to irradiation losses 
dominating surface oscillation behavior in the considered range of nanosphere dimension. Next we apply such a model 
to explanation of experimentally observed giant increase in photo-effect efficiency due to near-field coupling of surface 
plasmons with electrons in substrate semiconductor, when large metallic nanospheres are deposited on the optically 
active photo-diode layer. This coupling creates very effective channel for energy transfer as it will be demonstrated 
below. 



B. Solution of RPA equations: volume and surface plasmon frequencies 

Eqs (|15p and (|16J) can be solved upon imposed boundary and symmetry conditions — cf. Appendix [X] Let us 
represent both parts of the electron fluctuation in the following manner: 

Spi(r, t) = n e [fi(r) + F(r, t)\ , forr< a, ^ 
5p 2 {r,t) = n e f 2 (r) + a(Cl,t)S(r + e - a), for r > a, (r — > a+), ^ ' 

and now let us choose the convenient initial conditions, F(r,t)\ t= o = 0, cr(fi, t)|t=o = 0, (£1 — spherical angle), 
moreover (1 + f\(r))\ r=a = f 2 (r)\ r=a (continuity condition), F(r, t)\ r=a = 0, j p(r, t)d 3 r = N e (neutrality condition). 
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We arrive thus with the explicit form of the solutions of Eqs (|151 fT6|) (as it is described in Appendix lA")) : 

(18) 



fl(r) = _fex|±l e -fer(a-r)l^£I ) for r < a? 

-k T (r-a) 



f 2 (r) = [k T a - ^|±I (i _ e -2 fcT a )] e_^p>^ for r > fl> 



where kx — \l 67r "" e2 = \ —r- For time-dependent parts of electron fluctuations we find: 
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(20) 



where j;(£) = -^h+i/iiO — the spherical Bessel function, Y; TO (f2) — the spherical function, ui n i — lu p ^J 1 + -gr^i — the 
frequencies of electron volume self-oscillations (volume plasmon frequencies), x n i — nodes of the Bessel function 
k n i — x n i/a, u>oi = ujptj 2j^pi — the frequencies of electron surface self-oscillations (surface plasmon frequencies). 
From the above it follows thus that local electron density (within RPA attitude) has the form: 

p(r,t)=p (r)+p neq {r,t), (21) 

where the RPA equilibrium electron distribution (correcting the uniform distribution n e ): 

f l _ f n e [l + /i (r)], for r< a, ( 



ie/j(r), /or r > a, r — > a+ 
and the nonequilibrium, of plasmon oscillation type: 

a (V - / n ^(r,t), /or r < a, . 

The function F(r,t) displays volume plasmon oscillations, while a(fl,t) describes the surface plasmon oscillations. 
Let us emphasize that in the formula for <j(fi, t), Eq. (|20p . the first term corresponds to surface self-oscillations, while 
the second term describes the surface oscillations induced by the volume plasmons. The frequencies of the surface 
self-oscillations are 



WM = < ^V^Ti' (24) 

which, for I = 1, agrees with the dipole type surface oscillations described originally by Mie 11 , dj i = 0J p /y/3, 

C. Surface plasmon frequencies for metallic nanosphere embedded in a dielectric medium, with e > 1 

In order to account for influence of dielectric surroundings on the surface plasmons in metallic nanosphere, let us 
assume that electrons on the surface (r = a+, i.e. r > a, r — > a) interact with Coulomb forces renormalized by the 
relative dielectric constant e > 1. Thus instead of Eq. (TT6]) one can consider its following modification: 

( 25 ) 
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(Eq. (fTS"|) remains in not changed form) . Solution of the above equation is of the similar form as that for the Eq. (fTB)) 
case, but with new surface plasmon frequencies: 



The frequency of surface electron self-oscillations, changed by the factor J i, can be reduced significantly in comparison 

to the vacuum case, as in many materials e is relatively big (e corresponds to its high-frequency limiting value, the 
same which is involved in a refraction coefficient). 

Our result for resonant surface plasmon frequency, ^p^/j^r+ry > does not reproduce, for dipole case 1 = 1, the 

classical Mie formul a 23 ' 34 , uj v , 1 Our frequency is lower than the Mie one, which corresponds well with the 
data indicated in Fig. 3 in Ref. 1231 . in which there arc presented resonance frequencies obtained within a more 
thorough (TDLDA) method and located also below corresponding classical Mie values for several dielectric constants 
of surrounding medium (e = 1 1.7, 1.95, 2.31, 3 for air, Ar, Kr, Xe, MgO, respectively). Our formula has the similar 
property as includes some quantum effects (RPA approach) in comparison to completely classical Mie derivation 3 ^. 



III. EVALUATION OF A DAMPING RATE FOR SURFACE PLASMONS 

The RPA semiclassical equations (|151 IT6"| for plasmon excitations reveal the form of oscillator-equation-type. Thus 
it is easy to include, in the phenomenological manner, attenuation of these excitations, via damping term t§j 8p g^' t ' ) 
which can be added to the left-hand-side of Eq. (fTS"|) and Eq. (fT6|) (assuming that the volume modes, i = 1, and the 
surface modes, i = 2, are damped with the attenuation times rW, respectively). Thus the time-dependent solution 
of such modified equation (TTSJ) attains the form as given by Eq. (fTO|) with the factor e _ */ r<1) , and shifted frequency 
ui' nl = ^J^ni ~ ( T (i))2 f° r the volume modes. Similarly for the surface plasmons [Eq. (|16[) ]. the attenuation leads to 

the factor e~*/ T<2> for the first part of the solution (f2"Tj|) (and simultaneously shifted frequency lu' 01 — y^p; — - ^^y )i 

while the second term of Eq. ([20| acquires an additional factor e"*/ 1 "' 1 ' (and shifted frequency uj' nl = \J^ni ~ ( T (i)) 2 )• 

There are various mechanisms of energy dissipation of plasmon oscillations in metallic nanospheres. Let us concen- 
trate on damping of surface plasmons described by -^y. The interaction with phonons, electrons and lattice defects 
contribute to the relaxation rate — with the —fey, which is determined by the mean free path of electrons in the 

nanosphere, reduced additionally in comparison to the bulk case by inelastic scatterings with the sphere surface. One 
can use the estimation 3 ^ —^j ~ j 1 - + ^ £L , where vp — Fermi velocity, Xb — an effective value of the mean free path, 

Tj B 

C — constant of order 1, a — nanosphere radius (for Ag, vf = 1-4 • 10 6 m/s, Xb — 57 nm, which for a = 25 nm gives 

-j2j = 8- 10 13 s _1 , while rather a femtosecond decay time agrees with the measurements on Ag nanoparticles 3 ^) . Note 
T i 

that decomposition of surface plasmons due to creation of particle-hole pairs (Landau damping) 19 ! 21 is efficient only 
for small clusters 2 ^. 

Another type of energy dissipation can be associated with the radiation decay. The far-field radiation (i.e. for 
distances much longer than the wave-length A>a) gives the contribution to the relaxation ~ 3 2 ^ e - i io\ ~ 1.6 • 10 s 

s _1 (for u>i ~ 5 • 10 15 s^ 1 ) per single electron due to the Lorentz friction^. If one multiplies it by the electron number 

3 UJ 2 

N e = -^-n e , n e — (in order to account for the probability of energy transfer from the total system), one can 

arrive at the value -^y — {^~r§~^ > which dominates over -^y for not too small sphere s 27 ! 28 ' 35 . This channel of 

plasmon energy dissipation allows for explanation of the surface plasmon oscillations behavior with growing a (as -^y 

scales as a 3 ) for nanospheres embedded in a dielectric medium (like in the water — as it is presented in the Tab. 1 for 
nanoparticles of gold) . The more precise derivation of the far-field radiation losses expressed by — ^y is presented in 

Appendix IB! leading to the same formula for -^y as given above. Note that if the attenuation rate ujit^ is closer 



to 1 then the attenuation induced shift of the self-frequency is greater, oj' = 0J\"Jl — (uit^)^ 2 and this behavior 
coincides with experimental observations^ 7 - (for — 3_ > \ tfig overdamped regime is attained without free plasmon 

oscillations). 
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Tab. 1. Comparison with experimental data£ for Au nanospheres in water 



nanosphere radius 


Qj 


50 nm 


40 nni 


25 n ni 


attenuation rate due to far-field radia- 
tion losses 


(2) 


1.51 


2.95 


12.09 


shifted self-frequency rate 




0.75 


0.94 


0.99 


red-shifted oscillation energy 


hoj' (theor.) 


2.16 cV 


2.70 eV 


2.87 eV 


red-shifted oscillation energy 


Tiuj' (exper.) 


2.16 cV 


2.26 eV 


2.36 eV 



For Au, Tiu) p = 8.57 eV, and surface plasmon energy %uj\ = 2.87 eV. This value of hui is estimated assuming a 
coincidence of theoretically predicted self-frequency shifted by attenuation with experimentally measured for a — 50 
nm; note that the discrepancy between the experimental red-shift and the theoretical one grows for smaller a, which 
is probably caused by the strengthening of an impact of — hr ~ 1/a at smaller a, resulting in decrease of the red-shift 

in comparison to its value caused by -^y ~ a 3 . This tendency at decreasing radius a seems to be confirmed also by 

measurements for silver clusters with small dimensions < 10 nm, which was reported in Refs 381391 

Note that for small metallic clusters, quantum spill-out of electron cloud beyond positive jellium causes a red-shift 
of resonance Mie frequency lowering, however, with radius growth (thus it is in fact a blue-shift with radius growth). 
The additional effect of polarization of ionic system (this effect is beyond jellium model) can lead oppositely to inverse 
frequency shift, though rather small. Some jellium oscillation corrections can be also accounted for as scattering with 
phonons and can be included in the effective time rate for damping via effective mean free path A b ■ The corresponding 
contributions are, however, significant rather for small system s 16 ! 18 (for summarizing of various effects caused red and 
blue shifts of resonance with radius growth, cf. also Ref. 40). We consider the range of sphere dimensions when 
radiation losses cause overwhelming contribution to damping and to the resulting red-shift of the surface plasmon 
resonance, cf. Fig. 1. In this figure it is presented comparison of damping contributions due to scattering effects, 
~ f^" + ~ an< ^ ^ ue t° radiation losses in dielectric surroundings, ~ a 3 . For a > 10 nm the latter channel clearly 
dominates. The radiation-caused red-shift grows strongly with the radius of the nanosphere similarly as it is observed 
in the experiment for range of sphere radii 25-50 nm (for Au)£. 

The next source of the attenuation of surface plasmons would be connected with the transport of dipole oscil- 
lation energy between nanoparticles due to the Forster-type coupling 36 - in the case of sufficiently dense location of 
metallic nanoparticles. Nevertheless, taking into account that for uniform nanoparticle distribution in the dielectric 
medium, the same energy rates simultaneously escape and arrive at particular nanosphere due to interaction with 
other nanospheres (nearest-neighbors), this coupling does not contribute to the relaxation time (at least for uniformly 
distributed metallic nanocomponents) . 

The situation changes, however, significantly if metallic nanoparticles are deposited on the surface of the semicon- 
ductor substrate. Then the near-field e-m energy transfer from oscillating dipoles (surface plasmons with I = 1) to 
the electrons in substrate semiconductor starts to be the dominant channel of surface plasmon dissipation. One can 
estimate the corresponding time-rate by the Fermi golden rule applied to the system of plasmons coupled in near-field 
zone with semiconductor substrate. One can consider two situations, the first one — with rapidly switched off external 
electric field, which excites surface plasmons, then gradually (with lowering amplitude of oscillations) transferring 
energy to the semiconductor, and the second one — a stationary state of plasmons (with constant amplitude) with me- 
diating role of plasmons transferring entire energy of incident photons to semiconductor. The latter case corresponds 
thus to a stationary solution of a driven and damped oscillator, while the former one to free damped oscillations. In 
both cases the damping rate is the same, as it corresponds to the same substrate in a near-field zone. For free damped 
oscillations the total initial oscillation energy (assessed in Appendix [B| is gradually lost with the time ratio -L?. It 

allows for calculation of the —751, which is presented in the Appendix O Utilizing the similar calculus as in Appendix 

|B]one can assess the value of the corresponding damping rate -757, assuming that the total energy loss of surface plas- 

mons is transferred to the semiconductor substrate with additional renormalization by a factor j3 lowering an efficiency 
of this channel {[3 is a phenomenological factor introduced in order to account for geometry-induced proximity type 
constraints imposed on the dipole near-field coupling of the nanosphere with underlying semiconductor layer). Thus 
it is sufficient to calculate the energy income in the semiconductor due to nanosphere near-field dipole coupling — as it 
is done in Appendix [C] (within the Fermi golden rule scheme). For this channel of surface plasmon energy dissipation 
we deal with the scaling of the resonance energy shift with the dot radius, similar as that for -^y, however, with 

possible correction induced by dependence on a (it may be important, as for the nanosphere located on the planar 
semiconductor surface one can use an estimation (3 ~ c^r ~ 10 -3 , [for a — 50 nm] where c is a constant, h is an 
effective range of near-field coupling) . The parameter f3 significantly grows in the case when the total nanosphere is 
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in the near-field contact with the substrate, i.e. when the nanosphere is completely embedded in the semiconductor 
medium. For nanospheres deposited on the real semiconductor surface, the parameter (3 is obtained through fitting 
the experimental data (cf. Tab. 2). 

Assuming stationary conditions (i.e., constant in time amplitude of the surface plasmon oscillations, which corre- 
sponds to a balance of the incoming energy of incident photons with the energy outgoing to semiconductor substrate) 
the relevant damping is governed by the near-field dipole interaction (for fl<A) expressed by the scalar potential 37 
with an amplitude Dq(u>), 

ip{R, t) = -^—,n ■ D (oj)sin{ujt). (27) 

The matrix element of near-field dipole interaction for the transition of a semiconductor electron from the state in 
the valence to the conduction band, assumed as ^ fj(f) (r, t) = (2ir)~ 3 / 2 exp \ik-r — iEi^{k)t/h] (i-initial, /-final, 
respectively) is calculated in Appendix [Cl (Eq. (|C5jl ). which leads to a probability of transition per time unit, 

dw = — ^ 3(4^3)2^2^ " (hu> — Eg), where Dq(u>) is the surface plasmon dipole oscillation amplitude, adjusted to the 
balance of energy income and outcome (via shift of the resonance for stationary driven and damped oscillations) . 
Taking into account that the number of incident photons in the volume V of semiconductor equals g^°^ and the 

volume rate of metallic components Co = N m (N m — the number of nanospheres), the probability that an energy 
of a single incident photon is transferred to the semiconductor via surface plasmons on metallic nano-admixtures can 
be expressed as (with Sw given by Eq. (|C5j) ): 

OA7 . feElVY 1 l3C e 2 Luf 2 (uj)AiTa 3 , „ 

q m = PN m Sw (j^-) = 1287T 5ft4 »V^p<(fiu - E g ), (28) 

where f(u>) = — , 1 is the amplitude of forced surface plasmon oscillations. 

In order to assess an efficiency of the near-field coupling channel one can estimate the ratio of probability of 
energy absorption in semiconductor via mediation of surface plasmons (per single photon incident on the metallic 
nanospheres) to the energy attenuation in semiconductor directly from a planar wave illumination (also per single 
photon). In the latter case the energy attenuation in the semiconductor per single incident photon is given by the 

formula for ordinary photo-effect, q — m*ueh s ^fi u> ~ E g ) 3 ^ 2 (cf. e.g., Ref. l4lh . The ratio Sat turns out to be of 

order of 10 5 (at a typical surface density of nanoparticles, n s ~ 10 8 /cm 2 ) which (including the phenomenological 

factor [3 and H — semiconductor photo-active layer depth) is sufficient to explain the scale of experimentally observed 
strong enhancement of absorption and emission rates. It should be noticed that -4y grows with (3 (cf. Eq. (|C9[) ) 

and would attain the critical value for overdamped oscillator (— ^ — = 1), which precludes surface plasmon free 

oscillations. 

Very high efficiency (even if diminished by (3) of the near-field energy transfer from surface plasmons to semiconduc- 
tor substrate is caused mainly by a contribution of all interband transitions, not restricted here to the direct (vertical) 
ones as for ordinary photo-effect, due to absence of the momentum conservation constraints for nanosystems — cf. 
Appendix [Cj The strengthening of the probability transition due to all indirect interband paths of excitations in 
semiconductor is probably responsible for observed experimentally strong enhancement of light absorption and emis- 
sion in diode systems mediated by surface plasmons in nanoparticle surface coverings^^ 7 -'^. 

In the balanced state of the system when the incoming energy of light is transferred to the semiconductor via 
near-field coupling, we deal with the stationary solution of driven and damped oscillator. The driving force is the 
electric field of the incident planar wave, the damping force is the near-field energy transfer described by the -^y 

(assuming that this dissipation channel is dominating). The resulting red-shifted resonance with simultaneously 
reduced amplitude allows for the accommodation to the balance of energy transfer to semiconductor with incident 
photon energy. The amplitude of resonant plasmon oscillations Dq(uj) is thus shaped by f(u>) = — , 1 

V /K-^)2 +4w 2 /(r (2) ) 2 



The extremum of red-shifted resonance is attained at w m = lo\\J \ — 2{ujit^) 2 with corresponding amplitude 

2 ) / ( 9. /, ,2 iv( 2 h-2 \ rpi- i -c. ■ _,_-,_„_4.:„ r , i f„ 1 It,. (A 2 )\2 



~ Tg / I 2y u)\ — (t\ ) 2 J. This shift is proportional to \/{lji{t^ ') ) and scales with nanosphere radius a similarly 

(diminishes with decreasing a) as in the experimental observations 7 - (note again that for 1/t[ 2 ^ the dependence on a 
is opposite [grows with decreasing a]). 
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In order to compare with the experiment let us estimate the photo-current in the case of metallically modified 
surface in relation to the ordinary photo-effect. The photo-current is given by /' = |e|iV(g + q m )A, where N is the 
number of incident photons, q and q m are probabilities of single photon attenuation in ordinary photo-effect^ and 



due to presence of metallic nanospheres, i.e., of q 



2x/2 



e y/ 2 



3tt 6 m*ujsh c 



(hu - E g f /2 (cf. Ref. EH) and q rn given by Eq. 

V 

; A — y~ + is the amplification factor (t^' the annihilation time of both sign carriers, t„( p ) the drive time 

for carriers [the time of traversing the distance between electrodes] ) . From the above formulae it follows that (here 
/ = I'(q m = 0), i.e., the photo-current without metallic modifications), 
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27.53- 10 6 c — 




(29) 



where c = 4 ™ an 



fljf, with n s the surface density of metallic nanospheres, H the semiconductor layer depth, 4>(x) 



(x 2 -i)i+4x 2 /x' 2 ~7 1 ' x = x i = r 3 ^Ij x g = Eg/ (Jy^x), hu>i — 2.72 eV, m n ^ is the effective mass of conduction 

band and valence band carriers (for Si, m* = 0.19(0.98) m and m* = 0.16(0.52) m, for light (heavy) carriers, band 
gap Eg — 1.14 eV, e = 12), m is the bare electron mass. 

The results are summarized in Tab. 2 and in Fig. 2, for various radii of the nanospheres, and reproduce well 
the experimental behavior reported in Refl7|. By x m we denote frequencies corresponding to maximum value of the 
photo-current (i.e., to maximum of I' /I). 

Tab. 2. Comparison with the experimental data- for Au nanospheres on Si layer 



a [nm] 


n s [10 8 /cm 2 ] 




Lo m = x m huji (theor) [eV] 


Tno m (exp) [eV] 


<f>(x m ) 




50 


0.8 


0.772 


2.09 


2.25 


0.84 


1.55 


40 


1.6 


0.951 


2.58 


2.48 


3.00 


1.9 


25 


6.6 


0.997 


2.71 


2.70 


49.42 


1.75 



(the best coincidence with the experimental data is attained at ft = 3.5 ■ 10 



50 z 



r) 



(a[nm]) 2 > 

In Fig. 2 an estimation of normalized photo-current, /'//, with respect to wave-length is presented, for three 
sizes of metallic nanospheres (Au) deposited on photo-active Si layer, with structure parameters as listed below (the 



proximity parameter (3 
Legend to Fig. 2 



3.5-10" 



;_50j_ 

(a[nm]) 



panel in Fig. 2 


radius a [nm] 


concentration n s [10 8 /cm 2 ] 


layer depth H [fim] 


left 


(A) 25, (B) 40, (C) 50 


(A) 6.6, (B) 1.6, (C) 0.8 


2 


central 


(A) 19, (B) 40, (C) 50 


(A) 6.6, (B) 1.6, (C) 0.8 


230 


right 


(A) 25, (B) 40, (C) 50 


(A) 1.5, (B) 1.5, (C) 1.5 


230 



10 



5 040 



makes possible a significant growth of 



As it was indicated above, the relatively high value of ^ - j.u jjj mr 
efficiency of the photo-energy transfer to semiconductor, mediated by surface plasmons in nanoparticles deposited on 
the active layer, by increasing (3 or reducing H (at constant n s ). However, because of the fact that an enhancement 
of f3 easily induces the overdamped regime — cf. Eq. (|C9j) . a more perspective would be thus lowering of H, the layer 
depth (cf. Fig. 2 (left), where a significant growth of photo-current with lowering of active layer depth H illustrates 
the surface character of the effect). The overall behavior of I'/I(ui) = 1 + q m /q calculated according to the relation 
(|29|) . and depicted in the central panel in Fig. 2, agrees quite well with the experimental observations presented in 
Fig. 4 of Ref. (cf. inset in the central panel of Fig. 2), both in position, height and shape of photo-current curves 
for distinct samples (the strongest enhancement is achieved for a — 40 nm at densities as indicated above, in the 
Legend to Fig. 2), though q m /q is probably overestimated as the q denominator would be greater for doped real 
semiconductor structure which was not taken into account in the present calculus, similarly as surface effects — all of 
these would change the q denominator as well as its energy dependence, especially for longer wavelengths, where the 
discrepancy between theoretical model and experimental data is noticeable. 



IV. COMMENTS AND CONCLUSIONS 



The presented analysis featuring semiclassical RPA-type approach to collective fluctuations in metallic nanosphere 
in 'jellium' model deals with two types of plasmons, surface and volume ones. Within this approximation the self- 
frequencies of surface plasmon modes are independent of the sphere radius (similarly as classical Mie frequency for 
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dipole surface oscillations). There are, however, also surface modes induced by the volume modes and frequencies of 
these volume-induced surface oscillations depend on the sphere radius, similarly as the self-frequencies of the volume 
plasmons (given by the dispersion relation u; 2 ; = w 2 (l + x 2 ll /(kTa 2 )), x n i — nodes of the Zth spherical Bessel function). 
The e-m response of the sphere consists of both resonance types, the surface and the volume ones. It should be, 
however, emphasized that exciting of the volume modes is limited by the nanoscale of the system resulting in almost 
uniform e-m wave fields for resonant wavelength (dipole approximation regime). The uniform over the sphere, dynamic 
electric field excites the surface plasmons but not the radius-dependent volume modes. Therefore one can conclude 
that the experimentally observed significant dependence of resonant e-m frequencies on the radius of nanoparticles 7 
should be addressed to more complicated phenomena than radius dependent volume modes. 

The shift of the resonance frequency (in particular of Mie dipole type oscillations) for small clusters (up to a ~ 2 
nm) was analyzed against various quantum effects in microscopic type approaches, mainly of TDLDA typ o 16 ' 18 i 20 
also within semiclassical approaches^ 7 -. All these investigations indicate a major component of the experimentally 
observed red-shift of Mie frequency due to quantum spill-out effect (via reducing of density of electrons, resulting 

in factor Jl - ^ for resonance frequency, where the spill-out volume AN, i.e., number of electrons outside the 

jellium edge, was the subject of various microscopic estimation a 16 ' ' 18 ' 20 ) . Described in that manner red-shift of the 
resonance turns out, however, insufficient in comparison to experimental dat a 17 i 18 ' 21 . It is lower than observed in the 
experiment for ultra-small cluster s 16 ^ . The contribution to red-shift, also important rather in small clusters, was 
obtained due to decay of plasmons for particle-hole pairs (Landau damping) 19 i 21 , which improved fitting with the 
experiment. Additionally it was predicted an opposite blue-shift due to multi-plasmon anharmonic contribution 24 . 
For larger cluster it was indicated that the dominating spill-out factor weakens as ~ 1 16,17,18,20 an( j f or considered 
in the present paper region of nanosphere dimensions (above 10 nm up to 60 nm) is negligible in comparison to 
predominant in this region of nanosphere radii radiation loss contribution (cf. also Fig 1). 

Note, that within microscopic approach TDLDA more difficult it is to explain microscopically the observed width 
of resonance peak than its position 1 —. The processes which would contribute are: spontaneous ionization, Landau 
damping (i.e., interference of collective state with particle- hole pair with similar energy), evaporation of ion, ion 
vibrations (above jellium model). For larger a, above 10 nm, the dominating contribution to peak width 4 ^ is caused 
by irradiation effects which are stronger than scattering contribution ~ + (and C — constant of order 1 and 
Xb — effective mean path, were estimated by many authors, cited in Ref. I lot in particular including Landau damping 
type fragmentation of collective excitation). 

The quantum spill-out, Landau damping and coupling to ion excitations (beyond jellium model), though important 
for small clusters, are thus rather weak for large nanospheres and contribute resonance shift for radius range above 
10 nm far lower than experimentally observed. Another possibility for radius dependence of resonant frequencies is 
connected with the interaction of surface plasmons with other components of the system, which leads to damping of 
these oscillations. A shift of the resonance for driven and damped oscillator depends on the attenuation rate, which 
scales with the nanosphere radius. We have analyzed various channels of surface plasmon damping. The most effective 
channel for the surface plasmon damping is connected with the dipole-type near-field coupling of the surface dipole 
plasmons with semiconductor substrate, on which metallic nanospheres would be deposited, e.g., in nanomodified 
diode-type systems. Due to nano scale of the spheres for this coupling the momentum is not conserved, which results 
in a strong enhancement of the interband transition probability (because all indirect electron transitions between 
valence and conductivity bands in substratum have to be accounted for, provided energy conservation alone). It 
agrees with the experimental data referred to a significant growth of the energy transfer from surface plasmons in 
metallic nanoparticles to substrate semiconductor. 

In order to include the damping of surface plasmons one can introduce a phenomenological damping factor r 
(attenuation time) to the oscillation semiclassical RPA equation for electron local density fluctuations. As the form 
of the resulting equation is of the damped oscillator type, thus attenuation causes the red shift in the resonance 

frequency, (J = \Ju 2 — ^5. For driven and damped stationary oscillations, the red-shift of a resonance takes place 

with maximal amplitude at ui m = \J u 2 ~~~~ rpi- This red-shift is dependent on the sphere radius, via radius-dependence 
of T. 

Energy transfer to semiconductor surroundings mediated by surface plasmons is so effective that it may easily 
cause overdamped regime for plasmon oscillations. This channel is, however, reduced (typically by three orders in 
magnitude) by proximity constraints. Nevertheless, for nanospheres deposited on the semiconductor surface even only 
small fraction of the near-field channel (~ ^3- ~ 10~ 3 , for a ~ 50 nm, h is an effective range of the near- field coupling) 
causes a strong damping of plasmons. If to embed the nanospheres in semiconductor medium the plasmon system 
would fall in the overdamped regime (uj\t > 1). 

In the case of a small contact of the metallic nanospheres with the semiconductor substrate or in the case of an 
absence of semiconductor surroundings, a significant contribution to plasmon attenuation is due to far-field radiation 
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and electron scattering effects. The radiation contributions to - scale with particle radius a as a 3 (for both far- and 
near-field channels, though in the latter case the proximity constraints, included in (3, would modify this dependence 
to the linear one), while for scattering contribution, - ~ -. Thus the total attenuation rate — ~ Aa 3 + Ba + C— 
(A, B, C constants). For relatively big spheres (a > 10 nm) the radiation channels prevail, while for smaller ones the 
scatterings would be also import an t 38 i 39 . 

The reported strengthening of photo- voltaic effects due to plasmonic concentrators (layer of metallic nanoparticles 
on active semiconductor surface with n s of order of 10 8 /cm 2 ), for instance: up to 20-fold increase of photo-current in 
Si with nanoparticles Ag (40 nm [2-fold increase], 66 nm [8-fold], 108 nm [20-fold]) 6 , indicate the significant role of 
near-field energy transfer growing with the sphere radius. Another observations confirm also the strengthening role 
of plasmonic oscillations for emission and absorption phenomena in semiconductor diode systems, e.g., 9- fold increase 
of emission from Si diode modified with nanoparticles Ag of elliptical shape 120x60 nm, and resonant emission shift 
after covering nanoparticles Ag with 30 nm layer of ZnS3*i£ and up to 14-fold increase of absorption with various 
metal nanoparticles: Ag (12 nm [3-fold], Au (10 nm [5-fold], Cu (10 nm [14-fold])^. Influence of dielectric coating 
is caused by the surface self-oscillation sensitivity to dielectric surroundings (as for metallic nanosphere embedded 

in a dielectric medium), luqi — lo p \J 2rpr^! an d f° r typical e ~ 10 it gives strong decrease of resonant frequencies by 

factor ~ 0.3. The best correspondence with the experiment is attained for reported strong dependence of extinction 
features with respect to nanoparticle size (located on the surface of Si), the shift of the resonant peak corresponding 
to the change of Au nanoparticles radius: 25-50 nm (stronger extremum for 40 nm) 7 and simultaneous enhancement 
of photo-current seem to be well described by our model. 

Some experimental data indicate, however, existence of competitive mechanisms. For instance, for active medium 
Ti02 the photo-current diminishes in a wide spectral region (excluding the UV range) for coverings with nanoparticles 
Ag (3-6 nm), while the same coverings on optically active organic medium (DSC — dye solar cell) lead to strong 
increase of the photo-current for 3 nm Ag, but to the decrease of photo-current for 6 nm Ag^. The competitive factors 
can be linked here with retardation of the carriers transfer, despite plasmonic strengthening, or with the destructive 
modification of a photo-sensitive substrate material, for too small nanoparticles (more convenient are probably greater 
nanoparticles 6 -) . 

Summarizing, in the presented model nanosphere surface plasmons couple with substrate charges (band electrons 
in a substrate semiconductor) via photon-less short range e-m dipole interaction with very quick timing (thus very 
effective) — as confirmed by time-resolved spectroscopy measurements 7 . The strong enhancement of the efficiency 
results from nanoscale- induced incommensurability leading to all indirect in momentum interband transitions, not 
allowed for interaction of band electrons with the original incident planar wave photons as in an ordinary photo-effect. 
The type of dipole coupling is connected here with a specific e-m field gauge in the vicinity of the nanosphere within 
the distance lower than the wave-length (thus 'inside' the single photon), crucially distinct than for the planar wave 
(in the latter case only vector potential can be used, which is impossible in the former case)^. The schematically 
described above scenario qualitatively fits with the experimentally observed behavior and elucidates the timing of the 
particular steps of the energy-transfer-processes including mediating role of metallic nanosphere surface plasmons. 
The relevant time rates can be estimated within standard quantum mechanical attitude of Fermi-golden-rule type. 
Thus the presented above RPA plasmon description supply the convenient and simple tool for further modeling and 
optimization of the metallically nanomodified solar cell structures, towards enhancement of their efficiency. 
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APPENDIX A: ANALYTICAL SOLUTION OF PLASMON EQUATIONS FOR THE NANOSPHERE 

Below we present a method of solution of Eqs (fl~5l [TB)) . For the solution of Eq. (fl~5)) we assume in accordance with 
the rotational symmetry: 8pi(r, t) = n e [fi(r) + F(r, t)] , for r < a (with an initial condition F(r, t) = 0|t=o). After 
substituting it into Eq. (fT5"]) we obtain: 

V 2 /i(r)- 4/1(0=0, 
The function fi(r) (nonsingular at r — 0) has thus the form: 

AW = ( e ~ kTr ~ ekTr ) - (A2) 
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where a-const., kx — y 6 ™° e ' 2 — y -^r (fcr-inverse Thomas- Fermi radius), uj p — y i7rr ^ e2 (bulk plasmon fre- 
quency). For F(r,t) we assume a single harmonics F(r,t) = F u (r)sin(uji), as for a linear differential equation and 
suitably to the initial condition. Thus from Eq. (|A1D we obtain: 

V 2 F u (r) + k 2 F u (r) = 0, (A3) 

cj 2 — o> 2 

with fc 2 = a ,„ p . A solution nonsingular at r = has the form: 

^(r)=Aj,(*r)y Im (n), (A4) 



where ^4-constaiit, = ^/7r/(2£)/ i+1 / 2 (£)-the spherical Bessel function [J/(£)-the Bessel function of the first order], 
Y; m (fi)-the spherical function (f2-the spherical angle). Owing to the quasiclassical boundary condition, F(r, t)\ r=a — 
0, one has to demand ji(ka) = 0, which leads to the discrete values of k = k n i = x n i/a, (where x n i, n = 1, 2, 3..., are 
nodes of ji), and next to the discretization of self- frequencies u>: 

^ = u l i 1 + j^) ■ (A5) 
The general solution for F(r : t) has thus the form 

oo Z oo 

F(r,t)=^2 X! X! ^WK fc n;r)Yi m (0)sin(w n ;i). (A6) 

l—Q ra— — l n—1 

The neutrality condition, J p(r, t)d 3 r — N e , with 5p2(r, t) = a(ui, t)5(a + e — r) +n e f2{r), ( e — > 0), can be rewritten 

a oo a 

as follows: — J drr 2 /i(r) = f drr 2 / 2 (r), J d 3 rF(r,t) = 0, f df2o"(f2,t) = 0. Taking into account also the continuity 

o a o 

condition on the spherical particle surface, 1 + fi(a) = /2(a), one can obtain: / 2 (r) = (3e~ kT ^ r ~ a ^ / (kxr) and it is 
possible to fit a (cf. Eq. (jA2|) and /3 constants: a = ^ = fc Tfl _ fcra+I (1 _ e - 2fcrQ ). From the condition 

a 

/ d 3 rF(r,t) = and Eq. (|A6j) it follows that Aoon = 0, (because of / d£lYi m (uj) = 4Tr8io8 m o). Therefore the 


internal-volume electron fluctuations in the sphere have the form: 
Pl (r, t) = 0(o - r) [n e + 6 Pl (r, t)} = 

(A7) 



n e S(a — r) 



00 i 00 

1 _ feTa+l e -Ma-r)_i_ (l _ e" 2 ^) + £ £ £ ^(i^)^!!!)™^^) 

/=! m— — l n—1 



Now let us solve the Eq. (fTB]) for surface electron fluctuations assuming <5p 2 (r,i) = n e / 2 ('") + + e 

a), /or r > a, (r — * a+). In order to remove the Dirac delta functions we integrate both sides of the Eq. |T 



with respect to the radius length (J r 2 dr...) and then we take a limit to the sphere surface, e — > 0. It results in the 



following equation: 

a2 a^) = lime £ ,_ o __2_ j dr a. I [3^ + £f (£7(0) t)s{a + e >^ r)+ n j 2 (r))] r 2 S(a + e - r)} 

I 

00 

<r(Q, t) J drr 2 6(a + e - r)£5(a + e' - r) - a 2 n e / 2 (a) 


2 00 00 . 

d r ^.2 r JO g(^l,*V(a+e -ri)+rt e / 2 (ri) 



2e F 
3m 



(A8) 



/ drr 2 S(a + e - r) & / dr x r 2 J d ^^^ 

a 

2 00 a 

-£ / drr 2 S(a + e-r)&J dnr? J d^n /^+h^) 


Note that the last term in the above equation describes a coupling between volume and surface plasmon excitations. 
The two first terms of the right-hand-side of the Eq. (|A8j) vanish in the limit e, e' — + 0, (e < e') due to the delta 

function properties, and the third term gives the contribution: p a " (3(1 + k^a). The fourth and fifth terms contribute 
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in the following manner: 



W ^„e =lm= "' oo / oo ^ (A9) 

T l=lm=-ln=l 



In the derivation of the fourth term we have used (for a < ri): 

9 1 d ^ a 1 — 

1=0 'i 



9a ^a 2 + r\ - 2ar 1 cos^ da j-* r l x ' ' '—. r 



E -j+r p i( cos "?) = E 4+rfl(co S7 ), (A10) 



i 

where Pi(cosj) is the Legendre polynomial [P/(cos7) = ^ffj E ^im(^)^m(^i)]' 7 ^ s an an gl e between vectors 

m=—l 

a = ar and ri, while in the derivation of the fifth term (a > ri): 

In the derivation of the last term in the Eq. (|A9[) we used the explicit forms of F(r,t) and fi(r), together with 

a oo 

the neutrality condition, J drr 2 fi(r) = — j drr 2 /2(f) = — + kxa). The above described procedure leads to the 

a T 

equation for the surface plasmons: 

2 / * 00 / 

Z=0m=-2 CAlOl 
00 ! 00 a 1+2 

+UJpfl e E E ^imn5|^i / im(^)/rf?'l^f+2i;(Ani? , l)sm(w n ii), 

tOm=-ln=l 

where Wq ; = uo 2 2fp[ ■ Taking into account the spherical symmetry, one can assume the solution of the above equation 
in the form: 

00 1 

a(Q,t) =J2 E Hm(t)Yi m (Q,). (A13) 

1=0 m=-l 

After substituting it into the Eq. (|A12[) we find: 



^1 = 0, for 1 = 0, 



u^qimit) + E Wp"- e -4imn^- / ^i+2 ji(k n iri)sin(u! n it) , /or Z > 1, 



n=l 



(A14) 



The solutions of the above pair of equations have the form: 
<7oo(t) = 0, (acc. to initial condition), 



(1 + 1)uj 2 a r l+2 (A15) 

qim{t) = -^sin(ui it) + A; mra j^rT^+it^ n e / dr 1 ^ T 2ji(k n ir 1 )sin(uj n it). 

n=l p nl 



And finally, 



oo / 



= £ E ilm(n)^«n(woit) 

Z=l m— — l 



(A16) 



E E E ^"'" MoP-f^z+iL^ ^t 51 )^ / dr i^j+23i(kniri)sin(uj n it). 

l = \m=-ln=\ p ' ni 



Additionally let us comment on the equation for the case of nanosphere embedded in the dielectric medium with 
e > 1 — cf. Eq. (|25ll . The solutions of this equations are the same as presented above, however, with the frequencies 



wo; modified as follows: ujqi = u p J 2I+1 



1 1 
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APPENDIX B: CALCULATION OF TIME RATE FOR FAR-FIELD DIPOLE-TYPE RADIATION OF 

NANOSPHERE SURFACE PLASMONS 



In order to estimate the attenuation coefficient due to far-field radiation losses one can consider damping of 
nanosphere plasmons rapidly excited by switching off the uniform electric field, E(t) — Eq[1 — Q(i)]. The corre- 
sponding oscillations of the local electron density can be described by the equations: 



d 2 6 Pl (r,t) 
dt 2 



2 dSpi(r,t) v F 2 r / ±\ 

^(Tj -Q t = -g-Afyi(r, t) - u p Spi(r, t), 



for r < a, and 



2 Sp 2 (r,t) 
dt 2 



+ = -£* V [f n e + Sp 2 (r, t)f] S(a + e- r) 



£fV fd s r l¥ ^ (0(a - ri )S Pl ( ri ,t) + §0(n - a)6p 2 (r u t)) + ^E r (t) S(a + e- r), 



(Bl) 



(B2) 



for r = a (e — ► 0). For E not dependent on r, the driving force E(t) enters to the second equation only and leads 
to the driven solution corresponding to the dipole surface plasmon oscillations 5 P 2 = 5 / io(^)?io(i)- Similarly one 
can conclude that for nanospheres the visible light does not excite nanosphere volume plasmons as within the dipole 
approximation the incident wave electric field is uniform all over the sphere, unless the dipole approximation does not 
hold (i.e., when a ~ A). 

For E(t) = Eo[l — O(i)] (the rapid switching off the constant electric field Eq) the solution of Eq. (|B2|) has the 
form: 



?io(t) = 



/ 47r en e 
3 muf 



7 En 



1, for t < 0, 



-t/4*\ for 4 > 0) 



(B3) 



where to^ 



l 

r (2) 



and oji 



h is undamped dipole self-frequency. 



It is easy to calculate the loss of the total energy of the system A = £(t = 0) — £ (t = oo), i.e., by taking into 
account both kinetic and potential energy of electron system. Only potential interaction energy of oscillating electrons 
contributes, and £(t) = const. + ^a 3 q 2 Q {t), [the time dependent part of energy is caused by interaction of excited 
electrons, J rf 3 r 1; d 3 r 2 MjllMa^MM^ , wit h £l , e 2 -> 0, e 1 > e 2 ]. For q 10 given Eq. (Bj one 

can find, 



A 



e" 3 4tt 



£(t = 0)-£(t = oo) = —a 6 — 



en e E 



(B4) 



since £{t) = const. + f^f a 3 (^f) 2 (cosuj[t + s ^^j 



(2) 



On the other hand, assuming that damping of oscillations is caused by far-field radiation, one can calculate the 
energy loss A using the Poynting vector II = j-E x B, with v = cj ' \fe. The scalar potential of the e-m wave 
emitted by the surface plasmon dipole oscillations: p(r,t) = eqi (t)Yi (fl)S(a + e — r) is of the retarded form, 



<t>{R,t) = 



P\r,t- ! 



d 3 r, and for R > a, <j>(R,t) 



8r> ^L " - , here 



e\R-r\ " '' " iJ ^ iui " " VV"i "J eRv'" dt 

f rp(r, t — ^)d 3 r. In our case of surface plasmon dipole oscillations 



_R 

R 



D t 



eq w t 



rY 10 (Q,)6(a 



r)d 3 r 



0,0, eq 10 t 



R 



and the dipole moment 



(B5) 



Rc dt ' 

ek~v I ' r )divj(r, t — ^)d 3 r = j^n J j(r, t — ^)d 3 r for the sphere and due to div(j(n ■ r)) = (n • r)divj + j ■ n, 



Similarly, for the retarded vector potential we find A(R,t) = j^—— — , [because of |j| = —divj(r,t), <fi(R,t) 



which gives dD ^$ t - - = / j {r ,t — d 3 r\. 

Hence, for far-field radiation of surface plasmon dipole oscillations we have 



B = rot A 



y/i A d 2 D 
c 2 R nX dt 2 ' 



(B6) 
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and 



1 dA 1 
E = — —-V<P = -=Bxn, 
c at y/e 



(B7) 



which corresponds to the planar wave, and II = 9t eV 3R2 j (© is the angle between £) and 11). Next, taking into 

account that ^ = §U - ds, one can find A = J jj£dt = ^ J ( ^t^ ) d*. For L> 2 = e^Tr/Sa^o as in Eq. 
(IB5|) . with gio given by Eq. ()B3|) . one can find in this manner, 



oo 

x / dt 



A-^-^- 



-1 



- — a \"™f) Wl) 



-i 2 



77W 



Sill U)±t 7— (5) — U 



(B8) 



-2t/r. 



(2) 



(2^ f2l f2 ^ 

The latter integral equals to T2 /4, which together with, (o^Tj ) 2 + 1 = (wir^ ) 2 , leads to the expression: 



.4 



6e?; 3 3 

Via a comparison with Eq. (|B4[) . we finally find, 



e 4?r 6 f en e^V\ 4 (2) 



(B9) 



(2) „ / V3C 
V1T2 I 



(BIO) 



APPENDIX C: CALCULATION OF DAMPING TIME RATE DUE TO NEAR-FIELD INTERACTION OF 
SURFACE PLASMONS WITH SEMICONDUCTOR SUBSTRATE 

For the near-field regime (A > R > a, A>a) the vector potential has the same form as previously for far-field 

dL)(t— — "1 

since only condition a ^> R was used for its derivation 3 - 7 -, A(R,t) = — g t — . In the near-field region the e-m 
field is not of planar wave type and both vector and scalar potentials are needed to describe it. The scalar potential 

attains the form 4>(R, t) — ~dw * p" ' , (due to the Lorentz gauge condition 37 , divA = — ). The resulting Fourier 

components of fields B^ and E u (i.e. for monochromatic D — Doe - '"^* - ")) can be thus represented in this case 



as2I: 



and 



e » = H D ° + £ - + * ( * • a) (-s - S + ») } «**• (c2) 

where we use the notation for the retarded argument, ioj (t — — ) = zwi — ikR. For near-field region kR <C 1 one 

can neglect terms with -jj and Assuming also that for near-field e lkR = 1 one can obtain thus B u — and 
Eu> = [3n(n • £>o) — £> ], which corresponds to dipole electric field. 
The dipole type near-field potential can be written as follows: 

<p(R, t) = —rn ■ D sin(ujt + a) = w + e iut + w~ e~™\ (C3) 
eR z 

where w + — (w - )* = j^2-^ eMn 1 -Do; one can confine Eq. (|C3[) only to the if" 1 " term corresponding to energy 
absorption in the semiconductor. Then, according to the Fermi golden rule, the transition probability per time 
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unit between states &iki(i",t) = (2tt) 3 ^ 2 exp [iki ■ r - iEi(k 1 )t/h] , ^ 2k2 (r,t) = (2tt) 3 / 2 exp [ik 2 ■ r - iE 2 (k 2 )t/h] 
(semiconductor electron states from the valence and conduction bands, respectively), equals 

w(k u k 2 ) =y\< fe il w+ l fc 2 > | 2 «(^i(fci) - E 2 {k 2 ) + hw), (C4) 

where < fci|w + |fc2 >= J^pjs I e§i eIQn ' — - ^ 2> — d 3 R. Taking z axis along the vector q = k 2 — k±, then 
q ■ R = qRcosQi, n ■ Dq — Dq(cosQcosQi + sinQsinQicosipi) (O is an angle between Dq and q). Hence, 

OO 7T 27T 

< ki\w+\k 2 >= (2 Ju 2is e ia D J dR J sinQ 1 dQ 1 J d^cosQcosQx + sinQsin@ 1 cos(j) 1 }e i < !Rcosei = r§$qe ia ^, [as 
6 o o o e 9 



Jcoseismeideie" cosei = ~i£2^], and probability of transition to(fci,fc 2 ) = (2 ^f°^ cos ' e ^(^i(fci) - E 2 (k 2 ) + 

Tiki). In order to include all possible initial and final states in semiconductor, the summation with respect to k\ and k 2 
has to be performed (including filling factors f{k\) ~ 1 and f{k 2 ) ~ 0, as well as absorption and emission of energy). 
In the result we arrive with the total transition probability in semiconductor per time unit 8w ~ J 4-r / 4%iy(fei, k 2 ) 
caused by dipole surface plasmon oscillations on the single nanosphere. 

Let us emphasize that due to absence of the momentum conservation for the near-field dipole coupling in a vicinity 
of the nanosphere all interband transitions contribute, not only direct ones as for the interaction with the planar 
wave. It results in strong enhancement of the transition probability for the near-field coupling in comparison to 
photon (planar waves) attenuation rate in a semiconductor in an ordinary photo-effect. 

For the simplest model band structure, E\(ki) — E 2 {k 2 ) + huj = x + y — 7, where x = 2m l , x — 2— r and 7 — huj — E g 

(Eg is the semiconductor band gap) the integration over wave vectors gives the formula for the total probability of 
the transition 



0P V P 1- Eg), (C5) 



where n = — /, " . 

Assuming now that the dipole plasmon oscillations correspond to the damped oscillations which were excited by 
the rapid switching off the uniform electric field (as in the Appendix [B| . E(t) = Eq(1 — B(t)), with the dipole- type 
solution for electron distribution given by Eq. (|F33|) . we have 



D(t) = [0, 0, Dae- 1 /^ cos(uj[t)Y w (Vl)5(a + e-r)] (C6) 

with 

D = ^E^a 3 : (C7) 
in comparison to the Eq. (|B3[) we have neglected here the second term '""("i*) f or t^u![ well greater than unity. 

W 1 T 3 

One can now estimate the total energy transfer to the semiconductor (assuming that the dominant channel of the 
dissipation is the near-field interaction with semiconductor substrate and neglecting here the small shift of w{ due to 
dissipation) 



A = P dwh^dt = = f- 3 V " \ \ 9 ±, (C8) 

where (3 accounts for the proximity constraints which reduce the near-field contact of the sphere with the semiconductor 
medium; for the case of nanospheres deposited on the semiconductor layer surface (3 ~ \ ~ 10~ 3 , for a ~ 50 nm (h is 
an effective range of the near-field coupling), for the nanospheres entirely embedded in semiconductor surroundings f3 
would enhance significantly. Comparing the value given by the formula (|C8[) with the energy loss given by Eq. (|B4|) 
one can find 
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For nanospheres of Au deposited on Si layer we obtain: 



(2) 
T 3 W l 



( — ) 

V H / 



= 0.0059/3 r v " * , (CIO) 



TO TO 



for light(heavy) carriers in Si, m n = 0.19(0.98) to, m p = 0.16(0.52) to, and E g = 1.14 cV, e = 12, Tjwi = 2.72 cV. 
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10 20 30 40 50 60 5 10 15 20 25 10 20 30 40 50 60 

Radius [nm] Radius [nm] Radius [nm] 



FIG. 1: Comparison of contributions to surface plasmon damping (upper curve) of scattering term, ~ + ^p, and (far-field) 

radiation losses- induced damping, ~ a 3 , Eq. ()B10|I . for large metallic nanospsheres (Au and Ag); for sphere radius larger than 
10 nm irradiation-induced-damping dominates (horizontal dashed line indicates 10 13 level) 
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FIG. 2: Normalized photo-current y(A) for various parameters^: (3 = 3.5 ■ 10 3 r^^jys , (left panel) H = 2 /im, a = 25 nm 

(A), 40 nm (B), 50 nm (C), with densities n s = 6.6 (A), 1.6 (B), 0.8 (C) xl0 8 /cm 2 , (central panel) H = 230 (j,m, a = 19 nm 
(A), 40 nm (B), 50 nm (C), with densities n s = 6.6 (A), 1.6 (B), 0.8 (C) x 10 s /cm 2 , (right panel) H = 230 fim, a = 25 nm (A), 
40 nm (B), 50 nm (C), with densities n s = 1.5 (A), 1.5 (B), 1.5 (C) xl0 8 /cm 2 , respectively; coincidence with the experimental 
data 7 is achieved in the central panel, inset reproduces the experimental data* 7 - 



